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This report consists of the attached manuscript titled, "Series 
Integration of the Diaphragm Cell Transport Equation When the Diffusion 
Coefficient Depends Upon Concentration", by Judith B. Cain and James K. 
Baird. The manuscript will appear as an article in the September 15, 1992 
edition of The Journal of Chemical Physics . The work reported in this 
article was submitted by Judity B. Cain to the Graduate School of the 
University of Alabama in Huntsville as part of an M.S. thesis in chemistry. 
The University awarded Ms. Cain the M.S. degree at its June 1992 
Commencement. 

Integration of the the diaphragm cell transport equation in series was 
required in order to measure the diffusion coefficient in the system 
succinonitrile/water. A report on the measurments, themselves, is being 
prepared for publication. 



Series Integration of the Diaphragm Cell Transport Equation 
When the Diffusion Coefficient is a Function of Concentration 


Judith B. Cain and James K. Baird 
Department of Chemistry 
University of Alabama in Huntsville 
Huntsville, AL 35899 


Abstract 

In the case of the diaphragm cell transport equation where the interdiffusion 
coefficient is a function of concentration, we have derived an integral of the form, 

t = Bo + Bl In (Ac) + Bi(Ac) + B 2 (Ac)\ ... 

where t is the time and Ac is the concentration difference across the frit. The coefficient, 
Bo, is a constant of integration, while the coefficients, Bl, Bi, B 2 , ..., depend in general 
upon the cell constant, the compartment volumes, and the interdiffusion coefficient and 
various of its concentration derivatives evaluated at the mean concentration for the cell. 
Explicit formulae for Bl, Bi, B 2 , ... are given. 


*Submitted in partial fulfillment of the requirements of the M.S. degree from the 
University of Alabama in Huntsville. 



1. Introduction 


The diaphragm cell is the most convenient device for measuring diffusion 
coefficients in liquids. 1 ' 4 The cell consists of two well-stirred compartments containing 
solutions of different composition on opposite sides of a membrane, which is usually a 
sintered glass disk. The stirring guarantees that the contents of each compartment are 
uniform, so that mixing is limited to transport through the frit. In operation, the cell is 
oriented so that the normal to the plane of the frit is parallel to the gravity vector. The 
heavier solution is placed in the lower compartment, while the lighter solution is placed in 
the upper compartment; this arrangement obviates gravitational convection within the frit. 
Because molecular diameters are small as compared with the size of the interstitial spaces 
between the sintered glass beads, the diffusive transport within the frit is the same as in 
bulk solution. 

In the absence of volume change on mixing, the center of volume within the frit is 
fixed with respect to the laboratory. Diffusion through the frit can then be referred to 
volume-fixed coordinates, in which reference frame a two component solution has but 
one diffusion coefficient. 5 

The diaphragm cell is a relative device, since it involves a cell constant which 
must be calibrated by analyzing data for a solution whose diffusion coefficient is already 

known. 1 ' 4 Once calibrated, however, the cell may be used to determine an unknown 
diffusion coefficient by following the time relaxation of the concentration difference 
across the frit. 

For the important case where the diffusion coefficient is a function of 
concentration, we have recently shown that the relation between the relaxation time, t, 
and the concentration difference, Ac, is of the general form 

t = B 0 + B L ln (Ac) + BriAc) +^cf + ... (1) 

The Bo-coefficient is a constant of integration, while the remaining B-coefficients depend 
in general upon the cell constant, the compartment volumes and the interdiffusion 
coefficient and various of its concentration derivatives evaluated at the mean 
concentration for the cell. ^ our previous work, however, we have obtained explicit 



expressions only for Bl, Bj, and B 2 . 7 Below, we extend our proof to derive a formula 
for the general B-coefficient. 

2. Solution of the Transport Equation 

For a planar system, Fick’s first law states 

J = - D (c)| (2) 

where J is the flux, c is the concentration, D(c) is the functional form of the interdiffusion 
coefficient, and x is the coordinate normal to the plane. Regardless of the initial solute 
distribution within the frit, diffusion through the frit rapidly achieves a steady state. 8.9 
that case, J is independent of coordinate, and Eq.(2) may be integrated across the frit to 
obtain 


J = . 1 I D(c) dc (3) 

J Cl 

In Eq.(3), ( is the effective thickness of the frit, and ci and C2 are the concentrations 
below and above the frit, respectively. 

We let Vi and V 2 be the respective compartment volumes. If diffusion proceeds 
from Vi into V 2 , conservation of mass specifies that 

V.^U-JA (4) 

and 


V 2 ^ = ja 

dt 


(5) 


where A is the cross-sectional area of the frit. Substitution of Eq.(3) into Eqs.(4) and (5) 
leads to the coupled equations describing the diaphragm cell: 


Vi 


J f c ' 

1 

C2 

J r c 1 

1 

co 


D(c) dc 


( 6 ) 


D(c) dc 


(7) 



Combined with the initial conditions, ci(0) and C2(0), Eqs.(6) and (7) specify an initial 
value problem satisfied by the two concentrations, ci(t) and C 2 (t). As a pair of coupled, 
first order ordinary differential equations, Eqs.(6) and (7) are solved by two integrals. 

2.1 First Integral 

If Eqs. (6) and (7) are added, we obtain 

V 'f + V ’f = ° <8) 


The solution to Eq.(8) is 

5 = Vici(t) + V 2 c 2 (t) = VicKOl + V^O) 
Vi + V 2 V!+V 2 


(9) 


where c, the volume-averaged concentration, is a constant of the motion. 


2.2 Second Integral 

We divide Eq.(6) by V i and Eq.(7) by V 2 and subtract to obtain 

= P f D(c)dc (10) 

Jc2 

where 

p = f(v; + vi) (U) 

is the cell constant. Eq.(10) is ordinarily considered as the transport equation for the 
diaphragm cell. 1-4 

To solve Eq.(10) for the difference, Ac(t) = ci(t) - C 2 (t), we introduce into the 
right hand side, c = c + y, dc = dy, and 


Eq.(10) then becomes 


w _ ViV 2 Ac 
Vi+V 2 



( 12 ) 


D (c + y ) dy 


(13) 



For electrolytes away from low concentration (in which regime the Debye- 
Huckel-Onsager theory requires terms involving ( c^ and for nonelectrolytes at any 
concentration, we can expand D( c + y) in a Taylor series about c : 


D(c+ y) = D(c) + X y k (14) 

k=l 

where D^(c) = (d k D^Vdc 11 )^ . Substitution of Eq.(14) into the right hand side of 
Eq.(13) followed by term-by-term integration gives the result. 


where 


djffc =- |3D(c)Ac 


oo 

1 + X 

k=l 


(15) 


_ D k >(c) (V^ 1 - (-Vi^ 1 
(k+l)!D(c) ( Vl+ v 2 ^ 

Note in the ease Vi = V 2 that bk = 0 for k odd. 


(16) 


The object is now to separate the variables, t and Ac, in Eq.(15) and integrate. If 
we introduce the dimensionless time 


x = |3D(c)t 

the reciprocal of Eq.(15) can be written 


(17) 


When 


dx 

din (Ac) 


oo 


+ 2 



08) 


ibkfAcf 

k=l 


< 1 


the right hand side of Eq.(18) can be expanded in a geometric series. 


The result is 


(19) 


dx 

dln(Ac) 


= 1 + Y (-l) m 

m=l 


00 

I bk(Ac f 

'k=l ) 


( 20 ) 



The multinominal theorem is stated by Abramowitz and Stegun^® in the form. 




Z (n;ai,a 2 ,...,a n )* (xi) a, {x 2 ) ai ...(x n ) ao 


In Eq.(21), 


( 21 ) 


(n;ai,a 2 ,...,a ! )* = n!/l a >ai!2 82 a 2 L.n«-a n ! (22) 

is a multinominal coefficient. The inner sum in Eq.(21) is carried out over all integral 
values of the aj ’s such that 


a i + ai +...+a n = m (23) 

ai + 2 a 2 +...+ na„ = n (24) 

The values of the aj s may be established by considering the various partitions of n; aj is 

the number of times that the integer, j, occurs in that partition, while m is the total number 
of integers in the partition. 


If we identify bk with xk/k, we may substitute Eq.(21) into Eq.(20). After 
interchange of summations with respect to n and m, the result is 


. dx = ^ 
Ac _ 



(-lfm! Z( n;ai,a 2 ,...,a n )*(b 1 )»*(2b 2 ^...(nb n ) a » 


(25) 


In eq.(25), the variables x and Ac are separated, and the right hand side can be integrated 
term-by-term. After reintroducing the variable, t, through Eq.(17), we obtain for the 
integrated form of Eq.(25) the result 


t = B 0 + Btln (Ac) + BnfAc)" 

n=l 

where Bq is the constant of integration, while 


b l = -(pd©)‘‘ 


(26) 


and 


(27) 



(28) 


Bn = (Bi/nn!) £ (-1JT m! £ (n;a,,a 2 , (2bj)»... (nbj"- 


m=l 


Eqs.(26)-(28) solve Eq.(10) and complete the derivation of Eq.(l). Eq.(27) defines Bl 
while Eq.(28) is a formula for generating the polynomial B-coefficient of arbitrary order. 

Abramowitz and Stegun list in their Table 24.2 all the partitions of n for n = 1 to 
10 and give values for the corresponding binomial coefficients, (n;ai,a2,...,a n )*. 10 On the 
basis of that table, it can be seen that the complexity of the double sum in Eq.(28) grows 
slowly with n; specifically, as n advances from 1 to 10, the number of terms appearing in 
Eq.(28) is 1, 2, 3, 5, 7, 11, 15, 22, 30, and 42, respectively. This integer sequence 
represents the number of partitions of n and has been tabulated by Abramowitz at 
Stegun 10 and by Sloane. 11 

3. Discussion 

Of physical interest is the determination of the values of D(c) and the 
instantaneous derivatives, D(^) (c). This can be accomplished by fitting t v^ Ac 
experimental data to Eq.(l) and treating the B -coefficients as least squares parameters. 
Once the B-coefficients have been evaluated, D(c) follows from Eq.(27) while the D(^^c) 
follow from Eqs.(16) and (28). Since any t vs. Ac data set is necessarily finite in extent, 
Eq.(l) must first be truncated so that the number of B-coefficients to be determined does 
not exceed the number of experimental data points available. Accurate truncation is 
always possible, however, since for small enough Ac, the sum of all terms dropped can be 
made less than the last term retained. 12 Having derived previously expressions for Bl, 
Bi, and B 2 , we have successfully carried out this procedure to determine D(c), DO)(c) 
and D( 2 )[c) for aqueous hydrochloric acid at c = 1 M. 7 

To the extent that the accuracy of the experimental data will allow, the higher 
derivatives can be evaluated using Eqs.(16) and (28). Specifically, by virtue of Eq.(16), 
each bfc is proportional to D^)(c); the b^ are in turn related to the B n by a triangularized 
system of linear, algebraic equations. Employing Eqs.(27) and (28), we have tabulated in 
Table 1 explicit expressions for the first four B-coefficients. This triangularized system 
can be readily solved for the bk in terms of the B n . The first four solutions have been 
listed in Table 2. 



When integrating Eq.(lO), we could have expanded D(c) in a Taylor series about 
any fixed concentration consistent with the solution properties of electrolytes and non- 
electrolytes. We chose c, however, because as is evident from Eq.(16) and Table 1, this 
origin makes B n< = 0 whenever Vi = V 2 and n is odd. This property may be exploited 
experimentally in situations where only the value of D (c) is desired; since in this case the 
number of polynomial coefficients is reduced by half, a finite t vs,. Ac data set may be 

fitted more accurately, because Eq.(l) can be extended to higher order in Ac than would 
be possible otherwise. *3,14 
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Table 1. Triangularized system of linear algebraic equations for the Bn-coefficients in 
Eq.(l) expressed in terms of the bk-coefficients defined by (16). The parameter X= (3D (c). 

a Bn 

1 bi/X 

2 (b 2 -bi 2 ) / 2X 

3 (b 3 -2bib 2 + br)/ 3A 

4 (b 4 - 2bib 3 -b 2 2 + 3bi 2 ^ - bi 4 ) / 4X 



Table 2. Equations in Table 1 solved for the bk-coefficients in terms of the B n - 
coefficients. The parameter X = PD(c). 

k bk 

1 XBi 

2 2X.B 2 -X 2 Bi z 

3 3XB 2 + 4\ 2 B 1 B 2 - 3X% 

4 4 XB 4 + 2 X 2 (3B 1 B 3 + 2B 2 3 - 2 X 2 Bi 2 B 2 - x\* 



